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Abstract 

<^ We discuss double-diffractive (double-elastic) production of the r/-meson in the pp — > pn'p re- 

' action within the formalism of unintegrated gluon distribution functions (UGDF). We estimate 

also the contribution of 7*7* — ► rj fusion. The distributions in the Feynman xp (or rapidity), 
transferred four-momenta squared between initial and final protons (ti, t^) and azimuthal angle 
difference between outgoing protons ($) are calculated. The deviations from the sin 2 ($) depen- 
dence predicted by one-step vector-vector-pseudoscalar coupling are quantified and discussed. The 
results are compared with the results of the WA102 collaboration at CERN. Most of the models 
of UGDF from the literature give too small cross section as compared to the WA102 data and 
predict angular distribution in relative azimuthal angle strongly asymmetric with respect to tt/2 
in disagreement with the WA102 data. This points to a different mechanism at the WA102 en- 
ergy. Predictions for RHIC, Tevatron and LHC are given. We find that the normalization, t\ i 
dependences as well as deviations from sin 2 ($) of double-diffractive double-elastic cross section are 
extremely sensitive to the choice of UGDF. Possible implications for UGDFs are discussed. 
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I. INTRODUCTION 



The search for Higgs boson is the primary task for the LHC collider being now constructed 
at CERN. Although the predicted cross section is not small it may not be easy to discover 
Higgs in inclusive reaction due to large background in each of the final channel considered. 
An alternative way 0,05 HI is to search for Higgs in exclusive or semi-exclusive reactions with 
large rapidity gaps. Although the cross section is not large, the ratio of the signal to more 
conventional background seems promising. Kaidalov, Khoze, Martin and Ryskin proposed 
to calculate diffractive double elastic 1 production of Higgs boson in terms of unintegrated 
gluon distributions 0,IE00]- It is n °t cl ear at present how reliable such calculations are. 
It would be very useful to use the formalism to a reaction which is easy to measure. Here 
we shall try to apply it to the production of 7/ meson which satisfies this criterium. 

Recently the exclusive production of rf meson in proton-proton collisions was intensively 
studied slightly above its production threshold at the COSY ring at KFA Jiilich and at 
Saclay 0. Here the dominant production mechanism is exchange of several mesons (so-called 
meson exchange currents) and reaction via S\\ resonance (Io| . 

In the present note we study the same exclusive channel but at much larger energies (W > 
10 GeV). Here diffractive mechanism is expected to be the dominant process. In Ref.fl^ 
the Regge-inspired pomeron-pomeron fusion was considered as the dominant mechanism of 
the rj' production. 

There is a long standing debate about the nature of the pomeron. The approximate 
sin 2 ($) ($ is the azimuthal angle between outgoing protons) dependence observed exper- 
imentally 0] was interpreted in Ref.fl^ as due to (vector pomeron)- (vector pomeron)- 
(pseudoscalar meson) coupling. To our knowledge no QCD-inspired calculation for diffrac- 
tive production of pseudoscalar mesons exists in the literature. 

Of course, one should worry about the origin of the hard scale which may justify the 
applicability of QCD perturbation theory. As soon as the mass of 77' is not sufficient for 
that, it can be reasonably large 77' transverse momentum (or tip) which would serve as a 
hard scale. Bearing this in mind, we will examine the QCD result in the whole kinematical 
region, which should be understood as a sort of continuation of perturbative result to the 
region where its applicability cannot be rigorously proven. 




Both protons survive the collision. 
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FIG. 1: The sketch of the bare QCD mechanism. The kinematical variables are shown in addition. 

In Fig^we show the QCD mechanism of diffractive double-elastic production of rj' meson, 
analogous to the mechanism of Higgs boson production. We shall show here that approxi- 
mate (~ sin 2 ($)) dependence is violated in the QCD- inspired model with gluon exchanges 
within the formalism of unintegrated gluon distribution functions (UGDF) and a distortion 
from this dependence can help to select the correct model of UGDF. For completeness, 
in this paper we shall include photon-photon fusion mechanism shown in FigO which was 
sometimes advocated as dominant mechanism at high energies. 




FIG. 2: The sketch of the photon-photon fusion mechanism. Form factors appearing in different 
vertices are shown explicitly. 



II. FORMALISM 



A. Diffractive QCD mechanism 




FIG. 3: Kinematics of exclusive double-diffractive r/'-meson production. 
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The kinematics of the process on the quark level is shown in Fig. EJ The decomposition 
of gluon momenta into longitudinal and transverse parts gives 

k = -x[p! + k ,t = X2P2 + k 0}t , h = xi(/c - Pi) + ki tt , k 2 = x 2 (p 2 + h) + k 2 ,t • (2.1) 

We take into account below that Making use of conservation laws we get 

ki + p[ =Pi~ k , k 2 +p 2 + k =p 2 . (2.2) 

Taking the transverse parts from these relations gives 

h,t = -{p'i, t + k ,t), k 2 , t = p' 2)t ~ k ,t ■ (2.3) 

Following the formalism for the diffractive double-elastic production of the Higgs boson 
developed by Kaidalov, Khoze, Martin and Ryskin 0, H, H, 0| (KKMR) we write the bare 
QCD amplitude for the process sketched in Fig.l as 2 

a a9" '9" ■ 2 f i2i t//( ; n \fg!l( Xl ^ x 'l^O,t^ki t t^l)fg^( x 2,^ 2 jko t ,k2 t ,t 2 ) 

M g pp l pv l = tn d ko,tV(h, fc 2 , P M )-* p 2 9 2 -(2.4) 

The normalization of this amplitude differs from the KKMR one 0, 0] by the factor i. 
The bare amplitude above is subjected to absorption corrections which depend on collision 
energy. We shall discuss this issue shortly when presenting our results. 

The vertex function V(ki, k 2 , P M ) in the expression (|2.4|) describes the coupling of two 
virtual gluons to the pseudoscalar meson. We take the gluon-gluon-pseudoscalar meson 
coupling in the form 

V a p{k\, k 2 , Pm) = Vjy Fg* g *^ v '(kl, k\) e^ va/ 3 fcf k 2 . (2.5) 

Normalization is such that Fg*g*^ v / (0,0) = 1. The normalization constant Vn can be ob- 
tained in terms of the partial decay width T(r}' — > gg) as 

v »-- K W%f- NL0 - (2 - 6) 

The same normalization was obtained for the QCD double-diffractive production of x mesons 
in 0. 

The gauge invariance requires kfV a p = k^Vap = 0. On the parton level using Feynman 
rules we have to "hook" the gluon-gluon-pseudoscalar meson coupling ()2.5|) to the quark 
line by contracting it with incoming quark momenta. Replacing the quark lines by the 
proton lines one can show that the vertex factor V(k\, k 2 , Pm) in the amplitude (J2.4j) has 
the following form 

V = (k - Pl ) a (p2 + k fV afS = —^V a p . (2.7) 

X\ x 2 

Using relations (j2.1)l . ()2.3j) and ()2.5j) this expression can be transformed to 

V = V N F g * g *^ v i(kl, k\) e^af, (pi + p 2 y (pa + k f{p' + k Q ,) a {p' t - k , t f . 



For a sketchy derivation of this formula starting from the parton level one may look to |16( . 
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In the c.m.s. system pi + p 2 = and (pi + ^2)0 = \/s. Since k 0;t _L [p' lt x p' 2t ] we have 

V = -V N ^fsF g * g *^ v ,(k{,kl)ei k i(jp[ t + kQ >t )i{tf 2t - k Q>t ) k p u . (2.8) 

Introducing a unit vector in the beam direction of ingoing protons in c.m.s. n\ = pi,j/|pi|, 
we get finally 

V = -V N S - F g * g ^ v ,(kl k 2 2 ) [(p' lft + k , t ) x (p' 2it - k 0it )] n . (2.9) 

Normalization of this vertex function differs from the KKMR one jH[ by the factor 
s Fg*g*^ n i /2. Form factor F g * g *^ v i can be relevant at some kinematical regions and it should 
be taken into account. Factor s/2 makes the normalization of the full bare QCD amplitude 
(|2.4jl consistent with canonical normalization of the cross section (see Eq. (|2.29j) ). 
Expression (J2.8|) can be also written as 

V(h,k 2 ,P M ) = V N S -F g * g «^,(kl,kl) ■ |k M ||k 2it | -sin(0), (2.10) 

where is the azimuthal angle between k 14 and k 2t . In our case, in contrast to vector-vector 
fusion to pseudoscalars, <f> 7^ $. This may cause a deviation from sin 2 (<E>) distribution. To 
better illustrate the deviation we write the gluon-gluon-pseudoscalar meson coupling in the 
following equivalent way 



V(k 1 ,k 2 ,P M ) = -VN2 F g*g*-+v'( k i,t> k lt) 
Ip'mIIp'^I sin($) - |k 0)t ||p' lit | sm(ip + $) + |k 0)t ||p' 2)t | sin(^) 



x 



(2.11) 



Here $ is explicitly the azimuthal angle between the transverse momenta of outgoing protons 
p' lt and p' 2t , if) is the azimuthal angle between k ,t and p' 2i (0 < if) < 2n). 
It is convenient to introduce the dimensionless parameters 

IM t |p'i, t | |p' 2 ,tl 
p = , s = , v - 



m v ' m v ' rrirf 



We take k^ t = — |k 0it | 2 , k\ t = — |k 14 | 2 , k 2t = — |k 2t | 2 and write differential d 2 k 0t as: 
d 2 k j = —\k Qtt \d\ko t t\di/j. Then we obtain finally 



s f du 
' — x 



2w d , [^sin($) - ^sm(j) + $) + ^sin(^)] 

[fi 2 + £ 2 + 2fj£ cos(ip + $)] [/i 2 + r/ 2 - 2/J77 cos(^)] ' 
x /g,i^( x i5 ^o,t> ^i,t> ti) f g£ { x 2i %2i k^ t, k 2 t , t2)F g * g *^ ri i(k 2 t , k\ t ) . 

We have to take into account that the dimensionless arguments in /°{ , /° 2 and Fg*g*^ v > 



h 2 k 2 
o,t 2 M 2 tf 2 o t ^ / , 



— 5- = — a —77 + 2 ur? cos (?/>) 
m 2 , 
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are the functions of integration variables fi and if). So now we clearly see there is no simple 
angular behavior like Ai ~ sin($). The angular behavior of matrix element is more com- 
plicated. Only in the limit &o,t ~~ * the sin($)-behavior with some modulated amplitude is 
restored. 

We can obtain some information about angular behavior of matrix element from prop- 
erties of the integral (|2.12j) . Obviously, we have periodicity of Ai in $ with the period 
27T. For $ = and $ = tt we have immediately Ai = 0. It follows from the oddness of 
the integrand for these values of $. This investigation is a good check of numerical results 
shown at Fig. (j7J). Of course, more detailed information can be obtained only after numerical 
integration of Q3 with concrete tactions f°» and F„^,. 

The objects fgf/(xi, x[, k^ t , k\ t , t\) and f°fj (x 2 , x' 2 , k^ t , k 2 t , t 2 ) appearing in formula (|2.4j) 
and ()2.12|) are skewed (or off-diagonal) unintegrated gluon distributions. They are non- 
diagonal both in x and k 2 space. Usual off-diagonal gluon distributions are non-diagonal 
only in x. In the limit 2:1,2 —* x' X2 , k^ t — > k\, 2t and tip — > they become usual UGDFs. 

Using the relations (|2.1|) and k\ — k 2 = Pm and s ^> |k 0; f| 2 , we obtain 

sx x x 2 = m%, + |p' M | 2 + |p' 2it | 2 + 2|p / 1)t ||p / 2it | cos($) = m\, + |P Afji | 2 . (2.13) 

The longitudinal momentum fractions are now calculated as: 

m 2 , + \P M ,t\ 2 s 

X\,2 = exp(±y), 

s 

x' 1)2 = x = M. (2.14) 

Above y is the rapidity of the produced meson. 

In the general case we do not know UGDFs very well. It seems reasonable, at least in 
the first approximation, to take 

f${xi,x^kl p klt,h) = y/fPVi,^)-f?\xi,%t)-Fi(ti), (2.15) 
f${x 2 ,x' 2 X,tAt^) = \J fi 2 \A^,t) ■ f?\^, t ) ■ Fi{h) , (2.16) 

where Fi(t\) and Fi(t 2 ) are usual Dirac isoscalar nucleon form factors and t\ and t 2 are 
total four-momentum transfers in the first and second proton line, respectively. The above 
prescription is a bit arbitrary, although it is inspired by the positivity constraints |3J] for 
collinear Generalized Parton Distributions. It provides, however, an interpolation between 
different x and k t values appearing kinematically. Our prescription is more symmetric in 
variables of the first and second exchange than the one used recently in [35| for Higgs boson 
production. 

The UGDFs above have a property that 

f(x,k 2 )^0, (2.17) 

if k 2 —>■ 0. The small-A; 2 region is of nonperturbative nature and is rather modelled than 
derived from pQCD. Usually the UGDFs in the literature are modelled to fulfill 



k 2 



-■ F{x,k;) -»• const (2.18) 

t 
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if fcf — > 0. It is sometimes more useful to use F(x, kf) instead of f(x, kf). 

When inspecting Eq. (|2.4|) and ()2.16|) it becomes clear that the cross section for elastic 
double-diffractive production of a meson (or Higgs boson) is much more sensitive to the 
choice of UGDFs than the inclusive cross sections. 

B. 7*7* fusion 

It was advocated in Ref. 0] that the pseudoscalar mesons production at small transverse 
momenta may be dominated by the virtual photon - virtual photon fusion. In the following 
we wish to investigate the competition of the diffractive mechanism discussed in the previous 
subsection and the 7*7* fusion mechanism. 




FIG. 4: Kinematics of exclusive 7*7* fusion mechanism of ?/-meson production. 



In this case averaged matrix element squared 

W^E^WWl 2 . (2-19) 
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In the most general case the Born amplitude reads: 



a 



MM 



9/J.fJ.' 1 -\ 2 77 / 1 1 \ 1 ai P 9 * 



{u{p' 2 , A / 2 )[F 1 (t 2 ) 7 l/ ± i— k iy ,F 2 {t 2 )]u(p 2 , A 2 )} . 



(2.20) 



Limiting to large energies (y/s ^> m v > + M N + M N ) and small transverse momenta t\ and t 2 
(\t\ <C 4M^) the matrix element for pp — > prj'p reaction via virtual photon - virtual photon 
fusion can be written as 

MHt;X « eF^t ^Y^ Tll^'ih, k 2 ) iP2+ t P ' 2Y - «*k(*a) , (2-21) 
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where 



^ (^15^2) — ^ e F-y*^*^^/ {ki, /c 2 ) • (2.22) 



In Eq. (j2.21|) -Fi(ii) and ^1(^2) are Dirac proton electromagnetic form factors. In the fol- 
lowing we have omitted the spin-flipping contributions related to the respective Pauli form 
factors. F^*^*^ v r is a respective electromagnetic off-shell form factor normalized to 

F 7 * 7 *^,(0,0) = i -^, (2.23) 

where f v i is the meson decay constant. Alternatively one can use the relation 

where only measured quantities enter. Inserting PDG values of experimental entries for 7/ 
we get |F 7 . r _^(0, 0)| 2 = 0.116 GeV~ 2 . 
Now we can write 

( Pl+p ' 1 )^;/^'(k 1 , k 2 )(p 2 +p' 2 y « -*e 2 ■ (2 S )F rr _ >v (fc 2 , A|) |k 1)t ||k 2 , t | sin($) .(2.25) 
Collecting all ingredients together we get 

l-M^'l 2 - 4, 2 e 8 ^^) |F 7 , 7 *^,(A; 2 , fc 2 )| 2 |k lit | 2 |k 2 , t | 2 sin 2 ($) . (2.26) 

The t\ and t 2 dependences of F 7 * 7 *_^^/ are the least known ingredients in the formula 
(|2.26|) . It is known experimentally only for one virtual photon. In the following we shall use 
two different forms of the form factors. The first one is inspired by the vector-dominance 
model: 

P (1.2 7,2 \ _ Fy*^*^ V >(0, 0) 

r7wl 11 2j ~ (l-fc 2 /m 2 )(l-A; 2 /m2)- ^ M > 



The second one is motivated by the leading twist pQCD analysis in Ref. [1| 

T? (1,2 k2\ _ -^7*7*->y(0 ; 0) (r> qo\ 

^f*r^n'^1^2) ~ (1 - (jfcj + fc2)/ m 2) • V-*) 



Both these forms describe the CLEO data |12( for one real and one virtual photon. 
C. Phase space and kinematics 

The cross section for the 3-body reaction pp — > prj'p can be written as 



da= — \M\ 2 -d 3 PS. (2.29) 
The three-body phase space volume element reads 

d3ps - dS^dg^S^ ■ w<» + » - * - ^ - ■ < 230 > 
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At high energies and small momentum transfers the phase space volume element can be 
written as 

d'PS = ^dhdtz^d&dQ 5 (s(l - £0(1 - 6) - m 2 v ,) , (2.31) 

where £i, £2 are longitudinal momentum fractions carried by outgoing protons with respect to 
their parent protons and the relative angle between outgoing protons $ G (0, 2tt). Changing 
variables (£1,^2) ~~ > (xp,m 2 ,) one gets 

d 3 PS = -^-dhdt 2 — d%F d$ . (2.32) 

2 H s^ + 4(m;, + |PM, t | 2 )/s 

It is more convenient for lower (but still high) energy to use variable Xp. However, at 
very high energies the cross section becomes too much peaked at xp ~ due to the jacobian 

J « - 1 - (2.33) 

Jx 2 F + 4m 2 ,/s Z "V 

and the use of rapidity y instead of xp is recommended. The phase space element in this 
case has the following simple form 

d 3 PS = -J—dt^dyd® . (2.34) 

2 8 7T 4 S 

If xp is used then 



«W» 1-^/4+^ =Ff- (2-35) 
In the other case when the meson rapidity is used then 

£ 1)2 «l-^exp(±y) . (2.36) 

Now the four-momentum transfers in both proton lines can be calculated as 

h , 2 = -fl-^M<. (,37) 

?1,2 4l,2 

Only if £ 12 = 1, t\ t 2 = ~P\/2f ^he l a tter approximate relation was often used in earlier 
works on diffractive production of particles. However, in practice £1^ 7^ and the more 
exact equation must be used. The range of t\ and t 2 is not unlimited as it is often assumed. 
One can read off from Eq. (j2.37J) a kinematical upper limit for which is: 

tifl < - (1 ~ 6 ' 2)2 m; . (2.38) 

In practice these phase space limits become active only for \xp\ > 0.2. The lower limits are 
energy dependent but are not active in practice. 

The Mandelstam variables for subsystems {proton^ + 7/} and {protonf, + rj'} can be 
expressed via other kinematical variables 

si >2 = £1,2) + 2ti i2 • (2.39) 

It is also checked if 51,2 > (m,/ + rn p ) 2 , but this limit is not active in the region of interest 
(central production). 
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III. RESULTS 



Before we present our results let us discuss some input parameters for our calculations. 
The partial decay width T(j]' — > gg) is not well known. Of course 

W -> gg) < Tj* » 0.2 MeV . (3.1) 

In the following we shall take the upper limit in order to estimate the cross section. 
The form factors responsible for off-diagonal effects are taken in the form 

Am 2 p ~ 2 - 79t i,2 

^{^1,2) = Ta — 2 7 — \7i + ' /nvi\2 ' (^-2) 

(4m 2 - ti, 2 )(l - ti,2/071) 2 

The k\ and k\ dependence of the form factor F a * a *^^ {k\ , is not well known as it is due 
to nonperturbative effects related to the internal structure of the r/-meson. In the following, 
in analogy to the 7*7* — > rf form factor, we take it in the factorized double monopole form 

Fg*g*-+v(ki, k\) = (1 - kf/A2 j {i _ kl/K2J . (3.3) 

We take k\ = —k\ t and k\ = —k\ t . The parameter A os ~ m p may be expected. In general, 
it can be treated as a free parameter in order to quantify the theoretical uncertainties. 

In the present work we shall use a few sets of unintegrated gluon distributions which aim 
at the description of phenomena where small gluon transverse momenta are involved. Some 
details concerning the distributions can be found in Ref. 24] . We shall follow the notation 



there. We wish to stress in this context that all the distributions considered give quite 
reasonable description of the HERA F 2 data. Only scale dependent distributions require a 
separate short discussion. 

The Gaussian UGDF is obtained as 

Fnaive(x, ^ , /4) = ^fi^" ( x , /40 " fGauss{l^)/^ , (3.4) 

Then the off-diagonal distribution is calculated as 

ff fA = JftVi, tf)mx u k 2 lt , M |) • F 1 {t 1 ) 



ff* = \/f^ Vim^, kf t , Ml) " F 1 (t l ) (3.5) 

The choice of the (factorization) scale here is not completely obvious. We shall try two 
choices: 

(a) i4 = m^ffi = m\, , 

(b) fij = Ql,f4 = m\, . i— i 

The first choice is similar as in j33j]. However, it is not obvious if the scale associated with 
the "hard" production (g*g* — > r]') can be used for the left part of the gluonic ladder where 
no obvious hard scale appears. Therefore we shall try also the second choice where we shall 
use Qq = 0.26 GeV 2 , i.e. the nonperturbative input for the QCD evolution in Ref. [42]. 

Let us start from the da/dxF distribution. In Fig. El we show the results of calculations 
obtained with several models of UGDF (for details see pij)- For comparison we show also 
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TABLE I: The measure of the skewedness S w / 2 °f azimuthal angle distributions for different 
UGDFs and different center-of-mass energies. In this calculation -0.5 < xf < 0.5, -1 GeV < t\p < 
0. 



W (GeV) 


S v/2 (KL) 


S n/2 (GBW) 


S w/2 (BFKL) 


29.1 


0.5990 


0.7889 


0.3615 


200 


0.5867 


0.6628 


0.3131 


500 


0.5629 


0.4983 


0.2990 


1960 


0.5019 


0.2622 


0.2814 


14000 


0.3870 


0.2283 


0.2617 



the contribution of the 7*7* fusion mechanism. The contribution of the last mechanism is 
much smaller than the contribution of the diffractive QCD mechanism. 

In Fig. El we present distribution in t\ and t 2 (identical) of the diffractive production and 
of the 7*7* mechanism (red dash-dotted curve). The distribution for the 7*7* fusion is much 
steeper than that for the diffractive production. 

In Fig. 0we show the distribution of the cross section as a function of the angle between 
the outgoing protons. In the first approximation it reminds sin 2 ($). A more detailed 
inspection shows, however, that the distribution is somewhat skewed with respect to sin 2 ($) 
dependence. This is due to the two reasons: 

(a) kinematical - caused by interrelations of integration variables due to finite phase-space 
limits (present also for the pomeron + pomeron — > rj' fusion model), 

(b) dynamical - caused by nonlocality due to the internal loop of the diagram shown in 
Fig. (the sin(4>) dependence is embedded only in the loop integration). 

In order to quantify the effect we define the parameter of the skewedness of the $ distribution 

as 

s " 2(w) = KW« ■ (3 ' 6) 

If we take more differential cross section in the definition above than S n / 2 = S n /2(W, y, t\, t 2 ) 
(or S n / 2 (W, xf, h, t 2 )). Of course -1 < S n / 2 < 1. For exact sin 2 ($) distribution S 7T / 2 = 
0. In Table 1 we show the skewedness parameter S 7T / 2 for our model for different initial 
energies W = 29.1, 200, 1960, 14000 GeV, relevant for WA102, RHIC, Tevatron and LHC, 
respectively. Generally, the larger energies the smaller S n / 2 . The phase space limitations 
cause only a very small skewedness. 

In Table 2 and in FigJS] we show energy dependence of the total (integrated over kine- 
matical variables) cross section for the exclusive reaction pp — > prj'p for different UGDFs. In 
the case of Gaussian UGDFs we show in Table 2 also results with the second choice of the 
factorization scale. The cross section with the second choice is much smaller than the cross 
section with the first choice. Quite different results are obtained for different UGDFs. This 
demonstrates once again the huge sensitivity to the choice of UGDF. The cross section with 
the Kharzeev-Levin type distribution (based on the idea of gluon saturation) gives the cross 
section which is small and almost idependent of beam energy. In contrast, the BFKL dis- 
tribution leads to strong energy dependence. The sensitivity to the transverse momenta of 
initial gluons can be seen by comparison of the two solid lines calculated with the Gaussian 
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TABLE II: Energy dependence of the cross section (in nb) for different UGDFs. The integration 
is over -4 < y < 4 and -1 GeV < t± 2 < 0. The second lines for Gaussian distributions are for the 
choice (b) of the factorization scale. No absorption corrections were included. 



UGDF 


29.1 


200 


1960 


KL 


0.2867(+0) 


0.7377(+0) 


0.4858(+0) 


GBW 


0.1106(+1) 


0.2331(+2) 


0.1034(+3) 


BFKL 


0.3279(-l) 


0.9205(+l) 


0.2188(+4) 


Gauss (0.2) 


0.6391(+3) 


0.1697(+5) 


0.2964(+6) 




0.3445(+l) 




0.2984(+3) 


Gauss (0.5) 


0.7389(+l) 


0.2705(+3) 


0.3793(+3) 




0.4199(-1) 




0.4094(+l) 


7*7* 


0.7764(-l) 


0.2260(+0) 


0.3095(+0) 



UGDF with different smearing parameter cr = 0.2 and 0.5 GeV. The contribution of the 
7*7* fusion mechanism (red dash-dotted line) is fairly small and only slowly energy depen- 
dent. While the QED contribution can be reliably calculated, the QCD contribution cannot 
be at present fully controlled. It is even not completely excluded that the QED contribution 
dominates over the QCD contribution in some energy window. 

At present it seems impossible to understand the dynamics of the exclusive 7/ production 
at high energy without a real measurement. The Tevatron apparatus gives such a possibility, 
at least in principle. In Fig|§]we present two-dimensional maps t\ x t 2 of the cross section 
for the QCD mechanism (KL UGDF) and the QED mechanism (Dirac terms only) for the 
Tevatron energy W = 1960 GeV. If \ti\, 1^1 > 0.5 GeV 2 the QED mechanism is clearly 
negligible. However, at |ti|,|t2| < 0.2 GeV 2 the QED mechanism may become equally 
important or even dominant. In addition, it may interfere with the QCD mechanism. 

In FigElwe show a two-dimensional map t x $, where t — t\ or t%. The bigger t, the 
larger skewedness with respect to $ = tt/2. The skewedness, which is almost independent 
of the beam energy, is a generic feature of the QCD mechanism, quite independent of the 
choice of UGDF. The observation of the skewedness seems to be a condition "sine qua non" 
for the confirmation of the QCD mechanism. 3 On the other hand, the observation of the 
sin 2 $ dependence, as for the lower-energy WA102 data, may be very difficult to understand 
microscopically. 

Summarizing, the reaction under consideration seems very promissing in better under- 
standing of the QCD dynamics in the nonperturbative region. 

In the case of Higgs (or heavy particle) production the large mass sets a hard scale. In 
the case of 7/ the mass is only about 1 GeV. In this case the hard scale can be obtained 
by selecting large transverse momenta or analogously large and |t 2 |. In principle, these 
variables could be controlled by measuring outgoing protons. As an example in Table 3 we 
have collected results for different windows in the (£1,^2) space for the KL UGDF. In this 
case the cross section is dropping down rather slowly with increasing \ti\ and |£ 2 |- However, 
this result depends strongly on UGDF used in the calculation. Therefore measuring the 



3 The absorption corrections should only increase the skewedness. 
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TABLE III: Cross section in nb for different cuts on t\ and £2 and -0.5 < xp < 0.5. The limits for 
the (— ti) and (—£2) windows in GeV 2 are given explicitly. This calculation was done for the KL 
UGDF. 





(0,1) 


(1,2) 


(2,3) 


(3,4) 


(0,1) 

(1,2) 
(2,3) 
(3,4) 


0.1587( 0) 
0.5437(-l) 
0.8436(-2) 
0.1558(-2) 


0.5437(-l) 
0.2257(-l) 
0.4033(-2) 
0.7978(-3) 


0.8436(-2) 
0.4033(-2) 
0.9033(-3) 
0.2089(-3) 


0.1558(-2) 
0.7978(-3) 
0.2089(-3) 
0.5753(-4) 



TABLE IV: Comparison of the cross section (in nb) for rf and r\ c production at Tevatron (W = 
1960 GeV) for different UGDFs. The integration is over -4 < y < 4 and -1 GeV < £i )2 < 0. 
The second lines for Gaussian distributions are for the choice (b) of the factorization scale. No 
absorption corrections were included. 



UGDF 


rj 


Vc 


KL 


0.4858(+0) 


0.7392(+0) 


GBW 


0.1034(+3) 


0.2039(+3) 


BFKL 


0.2188(+4) 


0.1618(+4) 


Gauss (0.2) 


0.2964(+6) 


0.3519(+8) 




0.2984(+3) 


0.2104(+4) 


Gauss (0.5) 


0.3793(+3) 


0.4417(+6) 




0.4094(+l) 


0.3008(+2) 


* * 

7 7 


0.3095(+0) 


0.4493(+0) 



cross section for different cuts on t 1 and t 2 would be a farther test of UGDFs. 

The formalism presented in the previous sections can be applied to a production of other 
pseudoscalar mesons. In Table 4 we have collected cross sections for r/ c meson integrated 
over broad range of kinematical variables specified in the table caption. Again we have taken 
an upper limit assuming r(?7 c — > gg) = r*°', which may be even more reliable in the case of 
i] c production. These cross sections are very similar to the cross section for rj' production 
and in some cases even bigger. The results with Gaussian distribution, a = 0.2 GeV and 
first choice of factorization scale seems excluded, as constituting too large fraction of the 
total cross section. This strongly suggests also that the analogous result for rf production 
must be questioned. This seems to open a problem of understanding the WA102 data in 
terms of the QCD mechanism discussed above. 

Our result shows that the measurement of double-diffractive double-elastic production of 
f) c should be possible. However, one should remember about very small branching fractions 
for different decay channels of r/ c . It is not clear to us at present if the missing mass 
technique could be used at the Fermilab Tevatron. This would help to avoid the small 
branching fraction problem. The results obtained with different UGDFs differ significantly 
Any measurement of the reaction would be then very interesting to estimate (or limit) 
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UGDFs in the nonperturbative region. 



IV. DISCUSSION AND CONCLUSIONS 

For the first time in the literature, we have calculated exclusive production of 7/ meson in 
high-energy pp — > pr]'p collisions within the formalism of unintegrated gluon distributions. 
This type of reaction exhibits an incredible sensitivity to the choice of UGDF, which makes 
precise predictions rather difficult. Measurements of this reaction, however, would help to 
limit or even pin down the UGDFs in the nonperturbative region of small gluon transverse 
momenta k t where these objects cannot be obtained as a solution of any perturbative evo- 
lution equation, but must be rather modelled. The usual procedure is to extrapolate the 
perturbative regime via a smooth parametrization. For most of inclusive reactions the de- 
tails of such a procedure are not essential. In contrast, for the reaction discussed here the 
extrapolation is crucial. 

The existing models of UGDFs predict cross section much smaller than the one obtained 
by the WA102 collaboration at the center-of-mass energy W = 29.1 GeV. This may sig- 
nal presence of subleading reggeons at the energy of the WA102 experiment or suggest a 
modificaction of UGDFs in the nonperturbative region of very small transverse momenta. 
Experiments on exclusive central production of 77' at RHIC, Tevatron and LHC would cer- 
tainly help in disentangling the problem. With some cuts on £1, t± and the reaction 
under consideration can be measured at the Fermilab Tevatron |44j . An exact evaluation of 
the experimental cross section for the Tevatron will be presented elsewhere. 

A reasonable description of the WA102 total cross section can be obtained with UGDF 
obtained by Gaussian smearing of collinear gluon distributions and rather small value of 
the smearing parameter o"o ~ 0.2 - 0.3 GeV, clearly pointing to a nonperturbative effect. 
If the parameter o§ is adjusted to the total cross section a reasonable description of the 
da/dxp around \xp\ < 0.2 is obtained simultaneously. This was not possible with the Regge- 
like (two-pomeron exchange) description of the reaction [12j which produced distribution 
too much peaked at xp ~ 0. However, our approach gives somewhat too much skewed 
(asymmetric around $ = vr/2) distribution in relative azimuthal angle between outgoing 
protons compared to the WA102 data. This model gives definite predictions at larger energies 
where the contribution of subleading reggeons should be negligible. Experimental data at 
different collision energies would verify the solution and shed more light on the dynamics of 
the 7/ meson production. At present the Tevatron apparatus could be used. 

Measurement at lower energies would be also interesting. Natural possibilities would 
be FAIR at GSI and J-PARC at Tokai. Such data could shed more light on the role of 
subleading reggeons which seems important in understanding the WA102 data. 

Due to a nonlocality of the loop integral our model leads to sizeable deviations from the 
sin 2 $ dependence (predicted in the models of one-step fusion of two vector objects). 

The 7*7* fusion gives the cross section of the order of a fraction of nb at the WA102 energy 
W = 29.1 GeV, i.e. much less than 1 % of the measured cross section. The 7*7* fusion 
may be of some importance only at extremely small four-momentum transfers squared. In 
addition it can interfere with the QCD mechanism, which is similar to the familiar Coulomb- 
nuclear interference for charged hadron elastic scattering. 

Finally we have presented results for exclusive double elastic r\ c production. Similar cross 
sections as for 7/ production were obtained. Also in this case the results depend strongly on 
the choice of UGDF. 
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In the present calculations we have calculated only so-called bare amplitudes which are 
subjected to absorption corrections. The absorption effects lead usually to an energy- 
dependent damping of the cross section for exclusive channels. At the energy of the WA102 
experiment W = 29.1 GeV the damping factor is expected to be of the order 5-10 and 
should increase with rising initial energy. A detailed analysis of absorption effects is left for 
a future separate study. 
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FIG. 5: da/dxF as a function of Feynman xf for W = 29.1 GeV and for different UGDFs. The 
7*7* fusion contribution is shown by the dash-dotted (red) line (second from the bottom). The 
experimental data of the WA102 collaboration 11J are shown for comparison. 
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FIG. 8: a tot a s a function of center of mass energy for different UGDFs. The 7*7* fusion con- 
tribution is shown by the dash-dotted (red) line. The world experimental data are shown for 
reference. 
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FIG. 9: Two-dimensional distribution in t\ x ti for the diffractive QCD mechanism (left panel), 
calculated with the KL UGDF, and the 7*7* fusion (right panel) at the Tevatron energy W = 1960 
GeV. 




FIG. 10: The cross section for pp — > prfp as a function of t {t\ or £2) and relative azimuthal angle 
<I> for the Tevatron energy W = 1960 GeV. In this calculation the KL UGDF was used. 
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